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SEMIEMPIRICAL METHOD FOR PREDICTING 


EFFECTS OF INCIDENT -REFLECTING SHOCKS 

ON THE TURBULENT BOUNDARY LAYER 

By S. Z. Pinckney 
Langley Research Center 


SUMMARY 


A method is presented for the calculation of the "boundary-layer velocity 
profile change and various thickness parameter changes across a two-dimensional 
incident -reflecting shock. The method employs a one -dimensional hypothetical 
flow, representative of the boundary layer, in conjunction with a postulated 
boundary -layer shock model. Curves, calculated by using the method, are pre- 
sented for boundary-layer velocity profile, boundary-layer thickness, momentum- 
thickness and displacement-thickness changes across an incident-reflecting 
shock for a Mach number range from 2 to 6, a Reynolds number range based on an 
equivalent length from 3*3 X 10^ to 3*5 X 1C>7, and Incident shock -turning angle 
range from 0° to 12°. Curves are also presented (empirically determined) for 
important incident -ref lecting shock — boundary-layer interaction lengths. An 
outline of a procedure which can be followed in the practical application of 
the method is also included. 

Comparison of the calculated velocity profile change (measured through the 
profile index change) and boundary-layer thickness change with experimental 
data gave good agreement. For example, the values of profile index change and 
boundary-layer thickness change are found to agree with data to within 10 per- 
cent or better for the Mach number range from 1.9 to 3-85, Reynolds number 
range based on an equivalent length of approximately 10*7 to 10®, and incident 
shock -turning angles up to 13 °. 


INTRODUCTION 


The optimum design and performance of inlets for hypersonic air-breathing 
propulsion systems present a number of problems of interest. The successful 
design of such inlets requires methods for closely estimating the growth of the 
turbulent boundary layer across shocks in the inlet compression flow field. The 
present paper presents such a method for two-dimensional, incident-reflecting 
shocks which is based on the usual conservation equations and employs certain 
empirical representations of parameters. In view of the empirical features, 
the method is regarded as a correlation procedure for determining the effects 
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of incident-reflecting shocks on the turbulent boundary layer as indicated by 
data available in the literature and by some unpublished data. 

A number of analytical investigations have been conducted on the changes 
in a turbulent boundary layer across incident-reflecting shocks which generally 
use one of two basic lines of approach. The first approach employs momentum 
integral methods of which the analyses of references 1 and 2 are examples. The 
second approach (for example, ref. 3 ) employs a one-dimensional hypothetical 
boundary layer representative of the real boundary layer together with a postu- 
lated boundary-layer shock model. Sufficient experimental investigations are 
available (see refs. 4, and 6 ) which provide measurements of changes in 
velocity profile and various thickness parameters through the shock-interaction 
region for evaluating the accuracy of the theories. 

In general, the analytical methods of references 1 to 3 predict the veloc- 
ity profile change accurately but do not predict satisfactorily the boundary- 
layer thickness change. The present correlation method is similar in approach 
to that of reference 3 and employs a one -dimensional hypothetical flow and a 
postulated boundary-layer shock model. The shock model of the present analysis 
permits the evaluation by empirical representations of the mass flow influx to 
the boundary layer due to turbulent mixing in the shock-interaction region, 
whereas the shock model of reference 3 does not include a means fqr the deter- 
mination of this mass influx. 

In the present correlation method the boundary layer upstream of the 
incident-reflecting shock is converted to an equivalent one -dimensional flow 
through use of the momentum, continuity, and energy relations by assuming 
exponential velocity profiles and the modified Crocco enthalpy distribution. 

A boundary -layer shock model is constructed employing conventional two- 
dimensional, weak shock theory along with empirical representations for impor- 
tant shock model dimensions which are based on schlieren and shadowgraph pic- 
tures and pressure survey data. Downstream of the shock model at the station 
where the theoretical pressure rise is realized in the real flow, the equiva- 
lent one -dimensional flow is converted back to the real boundary-layer flow 
and yields the velocity profile index and the velocity, momentum, and displace- 
ment thicknesses at this station. The empirical representations are required 
because turbulent boundary layers are not amenable to theoretical evaluation. 

The correlation method covers the ranges of independent variables avail- 
able in the experimental data as follows: adiabatic -wall upstream Mach numbers 

from 1.0 to 6 . 0 , incident shock-turning angles from 0 ° to 12 °, upstream values 
of velocity profile index from k . 6 to 11. 5, and upstream Reynolds numbers 
(based on an equivalent length) from 3 A X 105 to 2 X 10®. In addition to the 
dependent variables previously mentioned, the correlation method provides the 
significant geometric lengths of the interaction region, and the paper also 
presents a means of evaluating the velocity profile development in the 
constant -pres sure mixing region downstream of the interaction region. 

A preliminary version of this paper was presented in reference J} the 
present paper presents further refinements of data correlations made possible 
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"by additional data and additional correlations of value to the solution of the 
overall incident -shock turbulent boundary-layer interaction problem. 


SYMBOLS 


Cp 

c p 

c v 

f(M e ) 

h 

H 

K 

^1* ^2 

L eq 

L u 

M 

m 

N 

n l 


average friction coefficient for 
specific heat of air at constant 
specific heat of air at constant 

effective Mach number function, 

static enthalpy, CpT 
boundary-layer shape factor, 5*/ 8 
mixing constant 

lengths identified in figures 8 and 9 


turbulent flow 

pressure 

volume 


Hell - 7 -=± M [2' 1/S 


i t P 


distance required for distorted velocity profile at station d to 
return to an equilibrium profile 

distance from inviscid intersection of incident shock with wall to 
most upstream point at which pressure rise is detectable 

distance from most upstream point the shock pressure rise is 

detectable on wall to point downstream of shock where pressure rise 
on wall is 0-99 of pressure ratio across shock 

Mach number 


boundary-layer mass flow per unit area, L 
free-stream mass flow per unit area, 

profile index value corresponding to an exponential velocity profile 
having same value of m/m^ as that of an actual profile 


exponent of an exponential velocity profile 


^ = /y\ 

v 7 V 




which fits 


1 z 

as much of an actual boundary- layer velocity profile as possible 
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N, 


9/<P* 


P 

Ap 


profile index value corresponding to an exponential velocity profile 
having the same value of cp/cp^ as that of an actual profile 

pressure 

theoretical overall static -pres sure rise across incident and 
reflecting shocks 


dynamic pressure. 


pV 2 


R 

*x 

*5 

% 

T 

V 

w 


Ax 

y 

Y i 

a 


„ t 

% 


P 

P T 


adiabatic -vail recovery factor for turbulent flow (assumed 
to be O.896) 

gas constant 

Reynolds number based on x 
Reynolds number based on 5 
Reynolds number based on 0 
temperature 
velocity 

f 5 

boundary -layer mass flow, / Pb^b 

J 0 


ay 


distance measured in direction of free-stream flow parallel to a 
plate 

increment in the x-direction 
perpendicular distance from vail 

distance from vail to intersection of vail shock and incident shock 
single shock-turning angle, deg 
extra compression angle, deg 
shock vave angle, deg 

shock vave angle for portion of incident and reflecting shocks 
betveen vail and intersection of incident and reflecting shocks 


ratio of specific heats (lA), C-, 


y c v 
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5 


v_ h 

thickness of boundary layer at point where = 1.0 

~L 

S-o, 

boundary-layer thickness as indicated in figures 8 and 9 

r s / p, V,\ 

displacement thickness, / 1 - ° p )dy 

J o V p i v V 

r 1 

boundary-layer energy parameter, / 

J 


5" 

&d 

5* 

n 
1 ii 
e 


*b 
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thickness of boundary layer at point where y— = 0.99 

l 


PbVb 3 ^(S 


free-stream energy parameter, 

p 5 p^Y / y 

momentum thickness. / — , — 1 - 

mass density 


total momentum flux per unit area in boundary layer 

a 


> p i( ] 


1 + 7Mj 


boundary -layer momentum parameter, / PbVbMf 

J n \o 




free -stream momentum parameter. 


Subscripts : 

A,B successive shocks 

b boundary layer 

d maximum pressure rise point on wall downstream of an incident- 

reflecting shock 

e equivalent 

eq equilibrium 

station where the boundary layer reaches an equilibrium velocity 
profile downstream of station d 


E 



i incompressible 

l local free stream 

t stagnation conditions 

1 station upstream of a shock — boundary-layer interaction where the 

pressure rise is first detectable on wall 

l/7 corresponding to an exponential velocity profile having as an 

exponent 1/7 

3 longitudinal station as identified in figure l(c) 

x axial location 

1 

Primed values are at station 5^ unless otherwise defined. 

ANALYSIS 


Analytical Boundary -Layer Shock Model 

If the supersonic flow over a plate were inviscid, the pressure rise on 
the plate generated by an intersecting incident -reflecting shock would be 
instantaneous (at station 2 of fig. 1(a)) and thus cover an infinitesimal dis- 
tance. In the real flow, however, a boundary layer exists next to the plate 
and contains low kinetic energy and low Mach number flow in its lower regions. 
The low kinetic energy flow does not contain sufficient momentum to traverse 
the pressure rise of the shock without gaining additional kinetic energy from 
the outer layers of the boundary layer. Thus, a spreading out of the shock 
pressure rise occurs to provide the distance needed to add the necessary 
kinetic energy by turbulent mixing. This phenomenon is illustrated by the 
pressure distribution of figure l(a) where the pressure rise produced by the 
shock starts at station 1 and a maximum value is reached at station d and it 
corresponds to the theoretical inviscid flow ! s shock pressure rise. As shown 
in figure l(a), there may be an inflection near the middle of the wall pressure 
distribution curve if the shock strength is sufficient to produce separation. 

The complicated shock pattern within the boundary layer produced by the 
spreading out of the shock pressure rise is illustrated in the bottom sketch 
of figure l(a) which is based on shadowgraph and schlieren observations; typi- 
cal shadowgraph and schlieren pictures obtained under contract NASl- 3^89 by 
General Dynamics Corporation are shown in figure l(b). The incident shock 
approaches the boundary-layer edge from the left and curves after entering the 
boundary layer to produce higher shock angles with higher turning angles at the 
lower Mach numbers of the boundary layer. The Increasing angle of the incident 
shock after entering the boundary layer is a result of the boundary-layer shock 
attempting to match the free-stream pressure rise. The thickening of the low 
kinetic energy regions of the boundary layer caused by the adverse pressure 
gradient results in a shock being initiated upstream in the vicinity of the 
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wall (wall shock) which turns the flow away from the wall. This shock crosses 
the incident shock, proceeds toward the free stream, and forms the reflecting 
shock. Near the edge of the "boundary layer, the wall shock is joined by addi- 
tional expansion and compression waves which originate within the boundary 
layer. In actual boundary-layer shock configurations, the extension of the 
free-stream incident shock to the wall (inviscid intersection with wall), sta- 
tion 2, is a convenient reference point for longitudinal dimensions. The 
length of the interaction region upstream from the reference point is and 

the total length Lt» 

Figure l(c) shows the selected analytical boundary-layer shock model to 
be similar to the one indicated by shadowgraph and schlieren pictures. In 
order to avoid an analysis involving curved shocks, the boundary-layer flow at 
station 1 of the analytical model is converted to one -dimensional equivalent 
flow through the use of integrated values of mass, momentum, and energy corre- 
sponding to the profile index N. The one -dimensional flow is then treated with 
the two-dimensional shock system of the analytical model by using existing two- 
dimensional shock tables. At an appropriate downstream station, the uniform 
equivalent flow is converted back to boundary-layer flow by the reverse process. 
Integrations involved in converting back and forth from boundary-layer flow to 
equivalent flow were made by assuming exponential velocity profiles and a modi- 
fied Crocco enthalpy distribution as presented in reference 8 . 

The orientation of the incident and wall shocks in the analytical model 
(fig. 1 (c)) is fixed by locating the intersection of the two shocks at a dis- 
tance from the wall Yq corresponding to an empirical correlation of data 
obtained from schlieren and shadowgraph pictures. This analytical shock model 
produces a void region (cross-hatched) adjacent to the wall which in the real 
flow corresponds to a low kinetic energy flow (in some cases separated). 

Up to station d 1 no mass flow has been added to the boundary layer; how- 
ever, data indicate that by the time the flow reaches station d, a considerable 
amount of mass -flow addition has occurred. In the attempt to determine a means 
for estimating the mass addition between stations d' and d, it was noted that 
if the mass-flow addition were assumed to be that required to fill the void 
region at station d T at conditions corresponding to the equivalent flow at sta- 
tion d ! , a total mass flow at station d* was obtained which experimental data 
showed to be 8 to 12 percent less than the measured mass flow at station d. To 
eliminate this discrepancy between calculated and actual mass addition, the 
strength of the wall shock was increased by stipulating a turning angle greater 
than the angle 05 by the small angle <xq. By increasing the turning of the 
wall shock by the amount cx^, the void region height was made to be a direct 
measure of the necessary mass addition. It was found that good agreement was 
obtained between experimental data and the computed thickness change when it is 

assumed that aq is equal to 2 ° for oc 1 ^ 3 ° and- equal to for a \ ~ 3 °- 

The assumption of a linear variation of cx^ for ^ 3° was made in order to 
keep the extra turning angle cc^ less than the free-stream incident shock 
turning angle a >i . 
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As a result of the preceding considerations, the wall shock of the analyt- 
ical model (fig. 1(c)) has a turning angle of a-^ + and the boundary-layer 

portion of the incident shock has a turning angle of ct^. The height of the 

void region at station d 1 is a direct measure of the mass-flow addition between 
station d T and station d. The equivalent flow at station d 1 is turned out from 
the wall by the amount 05. The expansion wave with turning a£ originates at 

the shoulder of the void region and produces flow parallel to the wall at all 
locations downstream of station 3 - Calculation of the turning angle 05 by 

the methods outlined in a later section also yields a value of the velocity 
profile index Nj (corresponding to equivalent flow at station 3)^ which 
experimental data indicated to be equal to N^. 

For simplicity in the present analysis all the mass-flow addition is 
assumed to occur in the mixing region downstream of station 3j as indicated in 
the sketch of figure l(c), although in the real case the mass addition starts 
in the vicinity of station d ! . Under this assumption the results of the analy- 
sis show that the values computed for oc^ by using the extra turning angle 

and the conditions at station 3 are correct and that is equal to N^. 

Therefore, in the region between stations 3 and d, the velocity profile 
remains constant while the mass-flow addition occurs. The boundary-layer 
thickness at station d 6^ is determined by incorporating the mass -flow addi- 
tion in a simple continuity relation between stations 1 and d. Other boundary- 
layer parameters of interest at station d are readily determined from and 

5 ^. In the calculation of 5 ^ and the effects of the wall friction were 

assumed to be negligible inasmuch as this condition seemed to be an unnecessary 
complication in view of the empirical features of the method. 

The assumption that the height of the void region at station d ! is a 
direct measure of the mass -flow addition is considered to be reasonable in view 
of the general character of the phenomena. The size of the void region is a 
direct function of the location of the intersection of the incident and wall 
shocks with respect to the wall; the location of the intersection in turn is 
determined by the effect on the boundary layer of the adverse pressure gradient 
induced by the incident shock. Furthermore, the height of the void region at 
station d ! relative to the boundary-layer thickness is a measure of the dis- 
tortion of the boundary layer of the real flow and thus a measure of the rate 
of turbulent mixing. 

The model static-pressure distribution for the equivalent boundary-layer 
flow is shown in f igure l(a) . As a two-dimensional shock model was assumed in 
the present analysis for a boundary -layer shock model, all pressure changes in 
the equivalent boundary-layer flow will be instantaneous and thus of infini- 
tesimal length. At station 1 the free-stream static pressure is converted to 
the static pressure of the equivalent boundary-layer flow which is slightly 
lower than the free-stream static pressure. The analytical shock model is then 
traversed by the equivalent boundary-layer flow until station d is reached and 
then the static pressure of the equivalent flow is converted back to the actual 
static pressure. The only criterion relative to static pressure that is incor- 
porated in the correlation method is that the overall static-pressure rise 
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■between stations 1 and d is the same as that for the real flow. The differ- 
ences, at stations 1 and d, between the magnitude of the static pressure of the 
actual boundary layer and the static pressure of the equivalent boundary layer 
are inherent in the method chosen for the conversion to one -dimensional flow. 


Definition of Characteristics of Undisturbed Boundary 
Layer Upstream of Shock-Pressure Rise 


The first step required by the present method is the determination of 
integrated values of mass, momentum, and energy corresponding to the profile 
index U of the undisturbed upstream boundary layer for the purpose of con- 
verting this boundary layer to equivalent one -dimensional flow. The preferred 
procedure for the determination of these integrated parameters would involve 
the use of a single simply determined parameter. The Reynolds number correla- 
tion of the profile index Up presents a method for evaluating such a param- 
eter; an example of this type of correlation is presented in figure 18 of ref- 
erence 9 . The velocity profile parameter Up is the exponent of an exponential 

velocity profile —■ = which is the best fit of the actual boundary- 


layer velocity profile. The correlation of the profile index Up in terms of 
Reynolds number R xe is presented in figure 2. The equivalent Reynolds num- 
ber R xe of each data point of figure 2 is based on the length of flat plate 
x e at local free-stream conditions required to obtain the momentum thickness 
of the station with turbulent boundary layer from the plate leading edge. The 
experimental data of figure 2 include data from references 9 to 1^ as well as 
some data obtained under contract UASl- 3^-89 with General Dynamics Corporation; 
these experimental data cover a Mach number range from 0 to 6. The curve 
shown in figure 2 was obtained from an empirical representation for incom- 
pressible equilibrium velocity profiles of reference 6 and seems to describe 
the compressible experimental data very well. In choosing experimental data 
for figure 2, data were chosen which were believed to be essentially for flat- 
plate undisturbed flow and thus also have equilibrium velocity profiles. Com- 
parison of all data reveals no Mach number effect on the profile index Up in 
the Mach number range considered. 


As shown in figure 2, experimental data exist from which the profile 
index Up can be obtained for a very broad range of Reynolds numbers. In the 
correlation method of the present analysis, however, the parameter U and not 
w L is used to determine the integrated values of boundary -layer mass flow, 
momentum, and energy in a manner similar to the method used in formulating the 
tables of reference 8. The profile index N corresponds to an exponential 


profile having a value of 







identical with that of the 


actual boundary layer. A curve representing the relationship between Up and 
U was obtained from an empirical correlation of the existing experimental data 
and is presented in figure 3 along with some example data points (obtained 
under Contract UASl- 3 ^ 89 ) • 
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Conversion to Equivalent One-Dimensional Boundary-Layer Flow 

An equivalent one-dimensional boundary-layer flow, for the purposes of the 
present analysis, is defined as a one-dimensional flow which has a value of 
total mass flow equal to the corresponding integrated values of the actual 
boundary layer; this equivalent one-dimensional flow has the same values of 
momentum and energy which correspond to a boundary layer with an exponential 
profile (exponent N) and the same mass flow as the actual boundary layer. In 
ensuing comparisons of the results from the present correlation method and 
experimental data, it was not found necessary to match the momentum and energy 
in the one-dimensional flow with that of integrated values from the actual 
boundary layer to obtain significant results. The expression for the equi val- 
ient one -dimensional Mach number derived in appendix A in a form similar to 
the mass momentum method of reference 15 is as follows: 


where 


’W ‘ 


u z 


jn 


1/2 


N 

TT~I 




_i 

% 


1/2 


1 + 






1 + 7M§ 



( 1 ) 


( 2 ) 


Equation (l) shows that the equivalent Mach number is a function of the veloc- 
ity profile index N and the corresponding integrated values of mass, momentum, 
and energy as defined in appendix A. In using equation (l) to convert the 
boundary layer upstream of the shock system (station 1 of fig. 1) to a one- 
dimensional equivalent flow, the value of N for equilibrium velocity profiles 
is obtained as described in the previous section. For nonequilibrium velocity 
profiles, the value of N must be determined by some other method. Values of 
the integrated boundary-layer parameters may be obtained from figures k to 6 
as a function of Mach number at the edge of the boundary layer M^ and the 
profile index N. The calculation of the curves of figures k to 6 was per- 
formed on a digital computer; the modified Crocco enthalpy distribution of ref- 
erence 8 was assumed in the computation of figures k to 6. Solutions for equa- 
tions (l) and (2) were also obtained in the computations and are presented in 
figure 7 in terms of Hi ~ as a function of N and Hi . 


Values of other equivalent flow parameters can be obtained by using M e , 
N, M 1, and other appropriate equations. For example, the equivalent static 
pressure p e can be obtained from the momentum equation (appendix A) as 
follows : 



1 + 7M£ 


JL 

*1 


(5) 
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Boundary-Layer Shock -Turning Angle and Profile Index Change 


After the upstream Boundary layer is converted to equivalent one- 
dimensional flow, the next step in the correlation method is to establish a 
procedure for computing the turning angle ce^ of the individual boundary -layer 
shocks. The criteria, from which the computation of is developed, have as 


their basis the matching of the overall static -pres sure rise through the bound- 

i / P ^3 N 

ary layer 


with the theoretical overall static-pressure rise through 

V^Wh 

the shocks and expansion outside of the boundary layer 


3 


At station 3 


n ly 


of figure l(c), the pressure changes have been completed in both the boundary 
layer and free stream and no mass -flow addition to the boundary layer has 
occurred in the shock model. The calculations for a-^ will be performed 

between stations 1 and 3 (fig* 1(c)) instead of between stations 1 and d. Sta- 
tion d is not a suitable downstream station as the mass flow added to the 
boundary layer in the mixing region between stations 3 and d cannot be deter- 
mined until a*b is known. 



The value of the overall static -pressure rise [ ■ - - — ) through the shocks 

Vl ’ 1 h 

and expansion in the free stream can be determined from two-dimensional shock 

,p U3\ 


tables (ref. 16 ), but the means of obtaining the overall pressure rise 




through the equivalent boundary layer is more complicated. The pressure changes 
in the boundary layer include pressure ratios associated with converting the 
real boundary layer to one -dimensional flow and back again. The overall pres- 
sure rise through the boundary layer is as follows: 



(*o 


As outlined in the previous section, the value of the ratio 


can be 


determined from the values of M ^ 2 . 8X1(1 %. are known. The net pres- 


sure rise through the boundary -layer shocks 


■e,3 

? e,l 


(■ 


which includes the wall 


shock with a turning of cc^ + %, the incident shock with a turning of 
and an expansion with turning of is a function of the unknown turning 


angle cx^. The pressure ratio 


cannot be evaluated until the station 3 
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equivalent Mach number M e ^j is known and M e ^j also is a function of ct^. A 
trial-and-error method of solution must be used where an arbitrary series of 


values of 05 are assumed. For each value of a^, the ratio 


■e ,5 


is deter- 


ge, 1 

mined from figures 7 and 5 and from equation (3). The value of is 

selected from the series of values assumed either by graphical or analytical 
procedures which satisfy the criteria 



The selected values of will correspond to specific values of M g ^ ^ 

and Nj where M g ^ is a result of the two-dimensional boundary-layer shock 

calculations and Nj is determined from figure 7 * Experimental data show that 
the value of Nj so determined is equal to N<q. 


Geometry of Analytical Shock Model 

The mass-flow addition to the boundary layer due to turbulent mixing 
through the shock system is essential in the accurate calculation of the 
boundary-layer thickness change across the shocks. The method chosen for the 
calculation of this mass-flow addition has certain empirical features; there- 
fore, a boundary-layer shock model as similar as possible to that of the actual 
flow is used. Schematics of this shock model are presented for a weak incident 
shock in figure 8 and for a strong incident shock in figure 9 * One the 

critical dimensions in the diagrams of figures 8 and 9 is the distance from the 
wall to the intersection of the incident and wall shocks; the intersection 
point is within the boundary layer in figure 8 and external to the boundary 
layer in figure 9. This dimension Yq, as well as other important dimensions 
of the shock model required to determine the boundary-layer shock geometry, is 
not amenable to determination by analysis. Therefore, an empirical representa- 
tion for Yq, as well as for two other important dimensions of interest 1^ 
and L-t, is developed by using data from references k, 5 > 10, 12, 17 to 22, 
data obtained at Ames and Langley Research Centers, and data obtained by 
General Dynamics under contract NAS 1 - 3^+89 • The dimensions 1 ^ and Lq are, 
as defined in figure 1, the upstream and total lengths of the pressure gradient 
imposed on the wall by the shocks. 


The determination of Y-i , or ^ 1 , is accomplished through the use of an 

81 

"sr 

empirically determined curve for = — and a family of empirically determined 

% 

Ln Y n 

curves for — . The correlation curve for — is presented in figure 10 
S 1 Lu 
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where 




is given as a function of 


Ap . 

*1,1' 


the parameter Ap is the theoreti- 


cal overall static -pres sure rise across the incident and reflecting shocks 
external to the boundary layer. No systematic trends with changes in Mach or 

Y 

Reynolds numbers are detectable for - — . The family of curves for the nondi- 

^ Lu 

mensional upstream length of the interaction region is taken from cross- 

b l 

plots of the aforementioned references and are presented in figure ll(a) . In 


the curves of figure ll(a). 


Ap 


is again the independent parameter and the 


curves are for given values of and Rq at station 1. To avoid confusion, 

no data points are shown in figure ll(a); the correspondence with experimental 
data is presented in figure ll(b). 


In figures 12(a) to 12(c) cross-plots of the nondimen sional length param- 
eter were obtained by using the data of references 4, 10, 12, 17 * l8, 

6 1 

19, and 22, and unpublished data obtained at Ames Research Center. The range 
of Reynolds number, Mach number, and turning angle of the experimental data 
used to determine the correlations of the shock-interaction length param- 


eter — is presented in table I. The points identified in figures 12(a) to 
£>1 

12(c) by symbols are not experimental points but points from cross-plotted 
experimental data. The limit curve shown for attached shocks in figure 12(a) 
was determined from figure 4 of reference 23 corresponding to wedges. Fig- 


ure 12(c) shows that extrapolation of experimental results obtained for 


l£ 

5 1 


to 


Reynolds numbers not covered by experimental data is not advisable. 


TABLE I. - RANGE OF VARIABLES USED IN EMPIRICAL 

DETERMINATION OF 4L CURVES 
5 1 


Reference 

Mach number 
range 

Turning angle 
range 

Reynolds number range, R$ 

Unpubli shed 
Ame s data 

2.58 to 4.02 

8° to 12.5° 

1.95 x 10 ^ to 7.5 X 10^ 

28 

2.95 to 4.80 

6° to 9 . 5° 

1.18 X 10^ to 1.60 X 105 

5 

2.97 

7° to 12° 

2.65 x io5 

4 

2-99 

4° to io° 

1.60 x io5 

15 

1.44 

1.5° to 6. 5° 

1.49 X 10 ^ 

9 

3.85 

4° to 10° 

5.13 x 105 

11 

2.99 

8.4° to 11.6° 

3.81 X 10 ^ 

22 

3.0 

6.6° to 8.1° 

3.14 x io5 

18 

3.0 

6.6° to 8.0° 

3.98 x io5 
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In summary, the determination of the important length parameters Yp, 1^, 
and Lp of the shock model is accomplished through the use of the empirical 
curves of figures 10, 11, and 12. The determination of these length param- 
eters, in conjunction with the equivalent one -dimensional flow properties at 
station 1 and the turning angles ap and provide sufficient information 

to determine shock diagrams similar to those of figures 8 and 9.? either analyt- 
ically or graphically, upon using conventional shock tables (ref. 1 6 ) to deter- 
mine the shock angles. 


Boundary -Layer Thickness Change 

The initial step required to determine the boundary-layer thickness change 
between stations 1 and d is the determination of the mass-flow addition to the 
boundary layer between stations 3 and d. (See fig. 1.) It was found that the 
height of the void region at station d 1 is a direct measure of this mass -flow 
addition to the boundary layer. Station d T is defined as the station at the 
point where the wall shock intersects the edge of the boundary layer or at the 
point where the incident shock intersects the edge of the void region, which- 
ever is further downstream. The mass-flow addition is computed by using the 
equivalent flow conditions at station d 1 in conjunction with the height of the 
void region at station d 1 . Therefore, the ratio of mass flow at station d to 
that at station 1 is given by 


w d 

p e,d' 

M e,d'| 

h + 2_l Mejd , 2 

W 1 


M e,l ' 



. 1/2 


(5) 


where all the values of the equivalent flow parameters are obtained by the 


method outlined in a previous section. The quantity 


5 d' 
5 1 


can be evaluated 


graphically or analytically by the methods outlined in appendix B. 


The boundary-layer thickness ratio ~ can be evaluated by the following 


equation: 
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where the values of all parameters were determined previously. Specifically, 
the values of 311(1 are ^^ermined from figure 4 corresponding, 

respectively, to Np and 3%. 

l4 


and 



Boundary-Layer Growth Downstream of Station d 


In the determination of the growth of the turbulent boundary layer through 
an incident-reflecting shock system, the development of the turbulent boundary 
layer should be considered from station 1 to some downstream station where the 
velocity profile has returned to approximately an equilibrium profile. The 
choice of this downstream station is necessary because at the point on the wall 
where the theoretical static -pres sure rise through the shocks has been realized, 
station d, the boundary-layer velocity profile may be very distorted and thus 
Na will be considerably lower than the corresponding equilibrium value. Thus, 
the next step in estimating the boundary-layer growth must include a means to 
determine the distance from station d to the station where the velocity profile 
has reached a near-equilibrium value. This step should also include a means of 
determining the velocity -profile development as well as the changes in impor- 
tant boundary-layer parameters through this mixing region. 


An empirical equation for the change in the boundary- layer shape param- 
eter in the mixing region downstream of station d is presented in reference 24 
for incompressible flow: 


H 

h i/7 


- l 




(7) 


It was assumed in reference 2k that the equilibrium profile would correspond to 
an N of 7 and that K has a value of 0.23 for incompressible pipe flow. 
Equation (7) presents a means for determining both the required length of the 
mixing region and the velocity -profile development through the mixing region. 
However, as shown in figures 2 and 3, the equilibrium value of N is a func- 
tion of Reynolds number and does not necessarily have a value of 7j therefore, 
equation ( 7 ) is modified as follows: 


H 



1 



-K~ 


(8) 


The equilibrium shape parameter (h ] corresponds to the equivalent Reynolds 

, x V 


number 


( R xe), 


n which can be evaluated from R Q n 

d S jd 


by using the procedure given 


in appendix C. The constant K has been evaluated from the data of refer- 
ences 5> 25 , 26 , 27 j 20, 24, and 28 and from data obtained under contract 
NAS1-3489; the results are presented in figure 13* The values of the shape 

parameter H ^which is equal to — ^ used in equation (8) to solve for the 

experimental values of K presented in figure 13 were obtained from figure l4 
by using the experimental values of N through the constant -pres sure mixing 
region. The equilibrium profile index used in equation (8) for the 

determination of experimental points of figure 13 corresponds to the Reynolds 
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number Rq of the most upstream experimentally investigated point of the 
constant-pressure mixing region. 

Equation (8) indicates that Ax may approach infinity before an equilib- 
rium profile is reached. In order to obtain a practical solution as to the 
necessary length of the const ant -pres sure mixing zone, it is arbitrarily 
assumed that the shape factor H is allowed to develop until it has a value 
equal to H e q + O.IO^H^ - H e ^^. By using this assumption equation (8) can 

be rearranged to give 


L eq = 2.302 . 

6 d K 

where L e q is a specific value of Ax, Equation (9) provides a means for 
estimating the length L e q downstream of station d required to reach a near- 
equilfbrium shape factor. The value of K corresponding to the value of N e q 
chosen for station d is determined from figure 13- The value of N e q chosen 
for station d is the value of the profile index which would he obtained on a 
flat plate with the same value of Rq as station dj the details of obtaining 
N e q and thus H e q are presented in appendix C. 

The longitudinal station at distance L e q downstream of station d is 
assumed to correspond to an equilibrium shape factor arid is designated sta- 
tion E. The increase in the boundary-layer thickness between stations d and 
E is produced by the mixing (in addition to the mixing between stations 3 and 
d) of a significant amount of free-stream flow with the boundary layer to pro- 
duce a near-equilibrium velocity profile at station E. This thickness change 
can be evaluated by assuming a constant -pres sure process and by omitting the 
effects of wall friction. Since with these assumptions there are no external 
forces on the boundary layer between stations d and E, there is no change in 
the integrated momentum of the boundary layer and thus no change in the momen- 
tum thickness 9. Therefore, the boundary-layer thickness calculation is as 
follows : 



The value corresponds to the equilibrium value of N at station d and 

its corresponding 


N, 


9/qy 


The value of equilibrium N at station d was pre- 


viously described in connection with equation (8) but the corresponding value 
of N cp/cp^ is obtained from the empirically determined curve of figure 15. 

(Example experimental data are also shown in fig. 15 •) The value of the pro- 
file index of figure 15 associated with each value of N and thus each 

actual velocity profile is chosen to be the value of N which corresponds to 
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JL = r 1 ^b! d /y\ 
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dj^-] identical with 


an exponential profile having a value of 

that of the actual boundary layer 
nonequilibrium zero -pres sure -gradient velocity profiles, the values of 

and 

of N and N m / rn along with the expression 


Since the curve of figure 15 also holds for 

/ 0 \ w 4 

^ J are then evaluated hy using figures 4 and 5 and- the correct values 


0 = m_ JP 

s cp z 


(11) 


As noted, the preceding calculation omits the growth in the momentum thick- 
ness 9 due to wall friction along the mixing length. If desired, this effect 
can be estimated by using standard methods for evaluating the friction coeffi- 
cient and the results are modified accordingly. 


The length estimated by equation ( 9 ) when added to the length of fig- 

ures 12(a) to 12(c) provides an estimate of the overall length of run of bound- 


D E 


ary layer disturbed by an incident-reflecting shock. The thickness ratio — 

6 d 5d 

evaluated by equation (lO), when multiplied by the value of jr- of equation (6), 


provides a means of determining the overall thickness change chargeable to an 
incident -reflecting shock. Also shown by this section is that two successive 
incident shocks, designated A and B, which might occur in an inlet compression 
field, for instance, should be separated by a distance D where 


D - (Lt " Lu) a + L eq + 


( 12 ) 


if an equilibrium velocity profile is to enter shock B. 


RESULTS AND DISCUSSION 


Profile Index Change 

The present method for calculating the change in velocity profile across an 
incident-reflecting shock is based on the assumption of exponential turbulent 
velocity profiles; the exponent N was selected as the profile index. The 
values of N used in the method correspond to a mathematically exact exponen- 
tial profile having the same — as the actual profile. Since turbulent 

boundary -layer velocity profiles are not exactly exponential, the accuracy of 
the results obtained by using this assumption is of interest. Figure 1 6 pre- 
sents a comparison of experimental values of profile index at station d of 
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references k, 5 , and 25 as well as the data obtained under contract ^.31-3^89 
with curves computed by the present method. A portion of the experimental data 
presented in figure l 6 has nonequilibrium upstream velocity profiles; these data 
are identified on the figure. The greatest discrepancy between data and curves 
is approximately 10 percent which is considered to be ample substantiation of 
the method. 


The changes in profile index parameters across incident -reflecting shocks 
have been computed by procedures outlined in previous sections for a systematic 
variation in values of the prime variables. These variables include incident 
shock-turning angle or pressure ratio, upstream Mach number, and upstream pro- 
file index Nq. The purpose of the computations is to illustrate in detail the 
effects of the independent parameters and to furnish values of parameters to be 
used in the application of the method. The computed values of the profile 
index % are presented in figures 17(a) to 17(c) for Nq values of 5, 7, and 
9, respectively. The profile index Nq curves of figures 17(a) to 17(c) are 


presented as a function of overall pressure ratio 


2 ,d 

? 2,1 


for curves of constant 


Mach number M^q; superimposed on figures 17(a) to 17(c) are curves of constant 
turning angle . Results show that the profile index when expressed as a 

function of pressure ratio is a strong function of Mach number. However, a 
study of the slopes of the shock-turning-angle curves shows that most of the 
Mach number dependence would be removed by expressing 'Nq as a function of ; 

P 2 ,d 


in fact, Nq was deliberately plotted against 


P 


2,1 


to spread the curves for 


constant Mach number for easy viewing. Results also show that the profile 
index Nq is a strong function of upstream profile index Nq; Nq increases 

p 


substantially as Nq increases from 5 to 9 for the same Mj q and 


2 ,d 


The computations for Nq do not extend below a value of 1.25 because 
according to reference 29, separation will always occur for lower values of N, 
at least for incompressible flow. An empirical curve for compressible turbulent 
boundary -layer separation under incident -reflecting shocks derived from refer- 
ence 10 is given in figure 17(b) for Nq of J. The method used for converting 
the data of reference 10 to the variables of figure 17 (b) is given in appendix D. 
This separation curve indicates that separated flow occurs somewhere in the 
shock system for values of Nq ranging from 2 to 5 whereas incompressible 
criteria of reference 29 indicate that the flow at station d would have 
reattached for these cases. 


Changes in Boundary -Layer Thickness Parameters 
From Station 1 to Station d 

The experimental boundary-layer thickness change data of references k, 5, 
17, and 21 as well as data obtained under contract NASI- 3^89 are compared in 
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figure 18 with curves computed by using the present method; figure 18 presents 
calculated and experimental thickness changes as a function of free-stream 
turning angle . Some of these data do not correspond to equilibrium 
boundary-layer velocity profiles at station 1 (fig- l); these nonequilibrium 
data are identified on the figure. In these cases the curve representing the 
present method was computed by using the actual values of Nq and Rg and 
satisfactory agreement was obtained. All other data of figure l8 correspond 
to equilibrium velocity profiles at station 1 . Agreement with data is within 
10 percent which is considered to be adequate substantiation of the method. 


The results of systematic computations of the boundary -layer thickness 
change across incident -reflecting shocks are given in figures 19(a) to 19(c), 

respectively, for Nq values of 5, 1 ? and 9 and R xe values of 3.3 X ICp, 


3.5 x lot and 3.3 X 10 ?; these curves are for equilibrium boundary -layer 
velocity profiles at station 1 . The calculations were limited to Mach num- 


bers M- 


■1,1 


up to 5*0 "because of the limitations of the empirical correlation 

5 - 


T I J 

of ~ of figure ll(a). In general, the curves of decrease initially, go 

5 i 5 i 

through a minimum, and then increase as ccq is further increased. The 


increasing values of ^ are caused by high rates of mass influx into the 

5 i 

boundary layer produced by the severe distortion of boundary-layer velocity 
profiles at high turning angles. The results show the .turning angle 


corresponding to minimum values of 



generally increases with Mach number 


M^ The results also show the effect of Reynolds number or N^_ on the 

$d 

boundary-layer thickness ratio — is small for equilibrium velocity profiles 

s i 


at station 1. 


A separation curve corresponding to that of figure 17(b) is shown on fig- 
ure 19(b); in general, separation occurs under the shock at higher values of 

bA 

turning angle than those corresponding to minimum values of — — . In inlet 

s i 

design,' however, shock-turning angles smaller than those corresponding to mint- 
ed 

mum values are generally of interest because of a need of higher total- 


pressure recoveries in the inviscid flow than could be obtained with the shock- 

turning angles corresponding to minimum ^ values. 

°1 


The results of computations of 


displacement thickness change 


1c 

5*1 


ness change curves of figures 19(a) 


the momentum thickness change 



and 


corresponding to the boundary-layer thick- 
to 19( c ) are presented, respectively, in 
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figures 20(a) to 20(c) and 2l(a) to 2l(c); these curves are also for equilib- 
rium profiles at station 1. The curves of figures 20 and 21 were calculated by 
using the following equations: 



( 13 ) 


The values of 


— are obtained from figures 19 (a) to 19 (c) 
6 1 


of — are determined from equation (ll) and the values of 
8 

from the following equation: 


whereas 


5 


are 


the values 
determined 


5* -| _m 

5 - "mj 


0A) 


The values of — used in equations (ll) and (1*0 are obtained from figure 4 
m Z 

by using, to determine station 1 values, the value of the profile index at 
station 1, Nq, and by using, to determine station d values, the value of the 

profile index at station d, Nq, obtained from figure 17- The values of -1 

n 

used in equation (ll) for stations 1 and d are obtained from figure 5 using the 
values of Nqj/q^ determined from figure 15 corresponding to these Nq and 

Nq values. 


0,3 ^ 

In general, the curves of and — — behave in a similar manner to 


0i 




8 d 


those of the thickness change -=■ curves; however, the minimum values of both 


9d 


5*. 


■£=■ and u occur at lower turning angles. This observation is important as 


9 


5* 


1 »"1 

shocks with turning angles smaller than those corresponding to the values 

obtained for the minimum values of — are of major importance in inlet 

5 1 

design. 
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Overall Development of a Boundary Layer 
Through an Incident -Reflecting Shock 


The intersection of an incident-reflecting shock with a turbulent boundary 
layer distorts the boundary-layer velocity profile and, as a result, the veloc- 
ity profile at station d is not an equilibrium profile; the value of is 


lower than the value of N e q* The equilibrium profile index chosen for sta- 
tion d, as already pointed out, corresponds to the value of profile index which 
would be obtained on a flat plate with the same value of Rq as station d. 
Appendix C, which presents the means of determining the chosen value of N e q 


for station d, shows that the resulting equilibrium profile index is a function 
of the free-stream stagnation temperature T^^. The chosen value of N e q at 

station d is used in the calculation of the overall changes in all the boundary- 
layer thickness parameters (S, 0, and 5*) from station 1 to the equilibrium 

station E downstream of station d. The chosen value of N e q at station d is 


also used in the calculation of the nondimensional distance 


^eq 

Si ' 


the distance 


required for the distorted velocity profile at station d to return to an equi- 
librium profile. Therefore, as some insight into the overall effect of the 
intersection of a shock with a boundary layer is of interest, computations were 


made for the overall thickness changes 

L, 


5 


eq 


0 


eq 


and 


8 * 
u eq 

Sp* 


as well as for 


the nondimensional distance 
310° K. 


J eq 

Sl 


Si 01 

for a free-stream stagnation temperature of 


Computations for — — were made by using equation ( 9 ) along with the 

5 1 

5rl 

values of — obtained from figure 19 ; the results of these computations are 
6 1 

presented in figures 22(a) to 22(c). The values of K (fig. 13) used in equa- 
tion ( 9 ) correspond to a stagnation temperature of 310° K. The results of the 

computations of show it behaves in the same manner as the thickness- 

change curves of figure 19 ; for example, the minimum values of each constant 

l»eq 

Mach number curve of — — occur at almost the exact turning angle of the 

5 1 

8p 

corresponding curve for 


A marked decrease in the value of 


J eq 


5 l’ ~ ™ * * Bp 

with increases in Mach number M^p or profile index Np ^or R X e) ( 


occurs 


Computations of the 
tion (l0) and the curves 


thickness ratio i were conducted by using equa- 
5 i 

8(3 

for — of figure 19 . These computations were 
5 i 
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conducted also for a free -stream stagnation temperature of 310° K and the 
results are presented in figures 23(a) to 23(c), As equation (l0) ignores any 
change in momentum thickness between stations d and E (equilibrium profile 


station), the values of the ratio 


eq 


6i 


0 d 

would be identical with those for — 


of figure 20. However, the Improvement in the profile from station d to sta- 
tion E does result in a change in the values of the displacement thickness; 
therefore, as a knowledge of the overall effect of the shocks on the displace- 


ment thickness is also of general interest, computations for 


eq 




were also 


>* 


made and are presented in figures 24(a) to 24(c). 
were made by using the following equation: 


The computations for 


eq 



(15) 


where the values of -g- are determined from equation (l4) and the values of 

are obtained from figure 23* The values of — used in equation (l4) are 
6 1 m l 

obtained from figure 4 by using, to determine the station 1 values, the value 
of the profile index at station 1 and by using, to determine the station E 
values, the value of the profile index N e q, incorporated in equation (l0). 


Upon comparing the results of the calculations for (fig* 23) with 


the results of the calculations for 


(fig, 19 ), it can be seen that a con- 


siderable increase in boundary -layer thickness occurs between stations d and E. 


This is not the case, however, when the results of the calculations of 


5 * 
5 * 


eq 


(fig. 24) are compared with the results of the calculations of 




of fig- 


ure 21; the improvement in profile between stations d and E causes a decrease 
in the displacement thickness. 


Practical Application Procedure for Changes in Profile Index 
and Thickness Parameters 

In order to utilize in a practical problem the various thickness param- 
eter change curves which are furnished, two assumptions have to be made as to 
the boundary-layer flow problem being considered. First, the boundary-layer 
flow has to be assumed to be turbulent; and second, the boundary-layer flow 
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upstream of the incident -reflecting shock interaction with the boundary layer 
has to have an equilibrium flat-plate (zero pressure gradient) flow velocity , 
profile. The calculation across the incident-reflecting shock when the fur- 
nished curves are used then follows these steps: 

(1) The momentum thickness at station 1 (fig* l) is calculated by the 
most accurate means available and the value of Rqi is determined. 

(2) By using the procedure outlined in appendix C, the Reynolds number Rgq 
is converted to R xe and the equilibrium profile index N eq ^q for station 1 

is determined from figures 2 and 3* The value of is then 

obtained from figure 15 corresponding to N e q q. 

( 3 ) The values of the boundary-layer thickness 5q and displacement thick- 
ness 5 q* at station 1 can then be computed through the use of 0 q along with 
equations (ll) and (l4) . 

(4) The change in the thickness parameters 5, 5*, and 0 across an 

incident -reflecting shock with an incident shock -turning angle is then 

determined from the appropriate set of curves (figs. 19 (a) to 19 (c) or 23 (a) 
to 23 (c); 20 (a) to 20 (c); 2 l(a) to 2 l(c) or 2*4- (a) to 24(c)) by the use of 
Neq^i; and q. The change in the profile index across an incident- 

reflecting shock of incident shock turning angle <x^ is determined from the 
set of profile index change curves (figs. 17 (a) to 17 (c)) by using N e( ^q, 

p ^ d * « 
pi,i 1,1 

( 5 ) The value of the thickness parameter of interest downstream of the 
shock, whether at station d or station E, is obtained through multiplication 

of its station 1 value times the change read from the appropriate set of curves. 

An example calculation using this step-by-step outline is presented in 
complete detail in appendix E. 


CONCLUDING REMARKS 


There are two important features of the present method of calculation of 
the changes in a boundary layer across an incident-reflecting shock which 
enhance its value over that of previous methods. The first of these features 
is that the profile index change, and thus the profile change calculation, 
requires no information or knowledge of the boundary -layer thickness change. 

The second of these features is that a boundary -layer shock model is assumed to 
be of a type consistent with the boundary-layer development and shock configura- 
tion revealed by shadowgraph and schlieren observations of the actual phenomena. 
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Utilization of the two features stated in the previous paragraph in the 
calculation of the profile index change and boundary-layer thickness change 
give good agreement with data. For example, values of the profile index change 
calculated at Mach numbers 1.9 to 3*0 and turning angles up to 11° are found to 
agree with the data to within 10 percent or letter; houndary-layer thickness- 
change data in the Mach number range 1.9 to 3-85, Reynolds number range (based 

on an equivalent length) of 10? to 10®, and for turning angles through the 
free-stream incident shock up to 13 ° also correlate to within 10 percent. 

The results of systematic computations of important boundary -layer param- 
eter changes from upstream of the shocks to two selected stations downstream of 
the shocks are presented in the form of curves which can be used for direct 
application. The two selected stations consist of the point where the theo- 
retical maximum pressure rise occurs on the wall and the point where the 
boundary -layer flow returns to an equilibrium profile. For these stations, 
systematic computations were made for the boundary -layer thickness, momentum 
thickness, and displacement thickness. In addition to the previously men- 
tioned computations, systematic computations for the profile index at the theo- 
retical maximum pressure point were made. Because of the limits imposed by 
experimental data for important empirical parameters, the present method of 
calculation is generally limited to a maximum Mach number of 5*0. 


Systematic computations were also made for the nondimens ional length 


parameter 


•^eq 

Sl * 


This parameter is defined as the length of plate required for 


the distorted boundary layer at the maximum pressure point downstream of the 
shocks to return to an equilibrium velocity profile. As was done in the 

L eq 

thickness -change calculations, the calculations of — — are also limited to a 

5 1 


maximum Mach number of 5*0* 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 20, 1965* 
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APPENDIX A 


DETERMINATION OF EQUIVALENT MACH NUMBER 


A one -dimensional equivalent boundary- layer flow Is derived as a function 
of the profile index N and free-stream Mach number M^ from equations for the 
total mass flow, momentum, and energy. The expression for the equivalent Mach 
number function of this equivalent one-dimensional boundary-layer flow as 
derived in the "mass momentum" method of reference 15 is 


where 


and 



(Al) 


(A2) 


(A3) 


In equations (A2) and (A3) the static pressure across the boundary layer is 
assumed to be constant and equal to the local free-stream static pressure. 

In reference 15 the total temperature was assumed to be constant across the 
boundary layer but in the present method the total temperature is allowed to 
vary across the boundary layer and, thus, the energy equation must be included 
in the analysis. The energy equation with the equivalent total temperature 
incorporated is 


6mT e,t c p = (“hJb 8 + \ 


where 


(A4) 


(mh)b 


m Z h Z 



Pb v b h b d /y\ 
PZ Vz h z \s) 
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and 

^ = nz ^ ( a6 ) 

After rearranging equation (A4) and assuming that the static enthalpy h is 

Vb /y\ l/N 

given by CpT and that — = , an expression for the equivalent total 

temperature is obtained: 


N 


N + 1 


L e,t 


= Ti 


M 2 JL 
1j 
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■hi 


(A7) 


Substitute equations (A2), (A3), and (A7) into equation (Al) and rearrange 
to obtain an expression for the equivalent Mach number function f(M e ) in terms 
of the profile index N and free-stream Mach number: 


M* 


f(M e ) = 


/m X 1 / 2 / W , 7 - 1 „ 2 n\ l/2 

■(=r) (p-T + — m z 


1 + /Mj 


2 JP_ 
<P l 


(A8) 


m cp ri 

The parameters — y and are defined^ respectively, as 


m 

m Z 


jp. 

fi 





PbV 

pi v i 3 




8 * 9 


8 8 


(A9) 

(A10) 

(All) 


The results of computations that were done for the parameters of equations (A9), 
(A10), and (All) are presented in figures 4 to 6, respectively. The calcula- 
tions of these parameters were performed on a digital computer by using the 
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V-h 

following expressions for — 

v z 

conditions : 


APPENDIX A 


and —2. given in reference 8 for adiabatic -wall 
p i 



The value of 7 is assumed to he 1.4 and the value of the adiabatic -wall recov- 
ery factor r is assumed to be O.896. 
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APPENDIX B 


EVALUATION OF BOUNDARY-LAYER THICKNESS 

5 1 

5 d 

The evaluation of the boundary-layer thickness ^ by using the boundary- 

layer shock models of figures 8 and 9 requires the values of the shock-wave 
angles and the corresponding turning angles of the boundary- layer shock models 

as well as the value of ^ from figures 10 and 11. The shock- wave angles and 

the corresponding turning angles of the boundary-layer shock models are pre- 
sented in the following table: 


Shock-wave angle 

Turning angle 


% 

1 — 1 

«\ 

-rQ 

ca 

% + a b 

h>,2 

% + «b 

K,2 

% 


a l + a b 

8,1 

a z 


6 d 

The calculation of — using the values of the shock- wave angles and the corre- 

6 i 

sponding turning angles of the boundary- layer shock models depends on which of 
the three cases is applicable. The three cases and the corresponding equation 
8' 

of — for each are: 

5 1 


Case 1 : For Y^_ < 8-^ where and Z 2 are indicated in figure 8(b), if 

v lo , 

In <= 2-,, then — is given by 


Z 2 _ Y l/ 5 1 , Sin ( P fr,l " ^b) 6 * 008 (Pb,2 - %) / Yj 


tan p. 


b,l 


sin 


^,1 6ln P b,2 


(Bl) 


and is given by 

5 i 
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Sp Sx pan p. 


V tan^ <2 - ( % + oi)] 


;ln (fb,2 ' ( a b + a b)] co3 ( a b + a b) 


Sln ®b,2 


Z 2 ) Zx> and Y-^ satisfy the conditions of case 1, the following expres- 


5' 

sion for is used: 

S 1 


&d = sin (< 3 b,l - a-b) sln (Pb,2 - a b) / Yi\ Yi 

8 1 sin p b ^x sin 3 b ^ 2 \ 5 1/ S 1 


Case 2 : For Z 2 < Zx and Yx ^ &x ^fig • 8j eqs . (Bl) and (B2) are used 


2 o 2 1 \ bq 

for — and — the following equation for — is used: 

5x b ± 5 -, 


5d , 3ln (Pb,i - °b) 3l °(%,2 - %) / _ 5i\ + 51 . ]1 - l 

Si sin p b ^x sin 3b, 2 \ 5 1/ 5 1 5 1 


tan cl! 


Case 5 : For Z 2 < Zq and Yq > 5q ( Z 2 and Zq are indicated in 

5 d 

fig. 9(b) ), the following equation for — is used: 

S‘L 

_ sln [ P Z,l - (°b + %)] sin [ P Z,2 - ( a Z + %)] _ Yx\ + Yl + l_l - l_2 + 


in sin[p^ 2 - (a z - %)] 


Bl/ Bl 


with — given by 


;ln [ P Z,l - K + %)] Sin [Pz,2 - ( a Z + %)] 
sin sin[p^ 2 - (a 2 - 

tan (c^ + c^) 


. IiU h - 


tan 8' .. 
b,l 
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1 -j 

and — given by 
5 i 


l ± s in [Pb, 2 " ( a b + %)] cos (% + ®b) 


sin P D ,2 


sin [p;,i - (®b + 

] sin Pz 

,2 - ( a Z 

+ 

sin sin If 

n 2 + ia l 

*Z,2 - ( 

a Z “ ®b) 


1^1 tan j 

P b,2 " 

(% + a b, 

r 


11 -& ♦& 


-> (BT) 
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APPENDIX C 


DETERMINATION OF EQUIVALENT FLAT-PLATE REYNOLDS NUMBER 


By definition, the flat-plate equivalent Reynolds number R xe for a given 
boundary-layer station is the Reynolds number based on the length of flat plate 
(zero pressure gradient) x e required for turbulent flow to produce the same 
R@, as the flow in question, at the local free-stream Mach number and total tem- 
perature. Based on this definition of R xe , and the following assumptions for 
the average incompressible friction coefficient and the ratio of compressible to 
incompressible turbulent friction coefficient, 

c _ . = ih , = 0.472 

[log 10 R 0i - log 10 -±- 



c f _ _e_ 

C Fj) i 9i 


the expression for R xe is given by 




Re 


c F c F } ^ 


'F,i 


(C2) 


(C3) 


The incompressible average turbulent friction coefficient C-p ^ is obtained 

from figure 25(a), which is based on equation (Cl), by using the value of Rq^ 
given by 


R 0i 


He 

Cp/Cp^i 


(c4) 


The ratio of compressible to incompressible friction coefficient - as 

u F,i 

obtained from reference 50 is given in figures 26 (a) to 26 (c) as a function of 
R@, Mj, and T^t • The flat-plate equilibrium flow profile index N ec ^ and 

equilibrium flow shape parameter (-^eq^ can ^hen determined by the use of 
figures 2, 5, and l4. 


4 

III 
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APPENDIX D 


DETERMINATION OF REYNOLDS NUMBER Rgn 


Reference 10 gives empirical curves, for boundary-layer separation, of 
the incident shock -turning angle as a function of Mach number and Reynolds 
number Rg»; the thickness 8 " corresponds to the point in the boundary layer 


where 


V-u 

^ = 0.99. 

z 


Thus, in order to utilize the empirical curves of refer- 


ence 10 to determine the separation curve of figure 17 (b), a value of Rgtt had 
to be determined corresponding to each Mach number and an Nq value of 7* The 
value of R 5 " is calculated by using figures 2, 3, 4, 5, 15 , 25(b), 27, and 
28(a) to 28(c) along with the profile index given in the text. Figure 27 shows 
the empirically determined relationship (along with some example experimental 

Oil 

data) between the equilibrium profile index N e q and the ratio fig- 


ures 28(a) to 28(c) give the ratio of compressible to incompressible friction 

Cp 

coefficient in terms of Rq^, Mj, and T^ ^ as obtained from 

°F, i 9 

reference 30. 


The value of Rq is determined from figures 2, 3; 25(b), and 28(a) to 
28 (c) along with the following relationships: 



Q 

The value of — is given by the relationship 
o 


= J3SL JP 

& - “z <Pj 


(Dl) 


(D2) 


(D3) 
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where the value of is determined from figure 4 by using the value of 

and N = 7 and the value of is determined from figure 7 ty using the 

value of along with the value of ^cp/cp-^ given hy the empirically deter- 
mined curve of figure 15- The value of is then given hy 


* 8 " 


^0 
0 5 

5 


(d4) 


where — is determined from equation (D3) and the value of 
5 

from figure 27 and corresponds to the value of N - 7- 


6 ^ 

5 


is obtained 


35 


1 

ii 



APPENDIX E 


EXAMPLE CALCULATION USING FURNISHED BOUNDARY-LAYER 
THICKNESS -PARAMETER -CHANGE CURVES 


The following example calculation which follows the step-by-step procedure 
as stated in the text and which utilizes the various thickness -parameter-change 
curves furnished is presented: 

(l) The boundary-layer conditions at station 1 are: 

M Z ,l = 3-0 

*1, t = 310° K 

9q = 8.265 X 10 meter 


P 7 t = 6.895 x io 5 ? - e ^ fcons 

l > t meter2 

R 9 I = 4.07 X 103 (as determined from ref. 23) 
Incident shock -turning angle a ^ = 8 ° 

( 2 ) Determination of R xe : Equation (C3) gives 


R. 


xe 


% 

C F C F,i 


= 3.3 x io fc 




For this example calculation, the value of 


c F,i 


is 0.66 and is 


determined from figure 26(a) by using M 7 =3*0 and 

, L ’ x 
R 01 = 4.07 X 10^. The value of Cp^ is determined from fig- 

tire 25 (a) by using equation (c 4 ) 



= 6.17 X 103 
CF 

c F,i 
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Determination of N eq ^p: The determination of N e( ^p is accom- . 

plished through two steps: (l) enter figure 2 with R xe = 3>3 X 10^ 
to obtain Np, p = 6.1; (2) enter figure 3 with Np^p = 6.1 to obtain 

^eq, 1 = T*0* 

Determination of (^ e q^ l) cpy cp : Eirfcer fig ure 15 with N e g^p = 7-0 

to obtain (N~ n i ^ = 7*1^-* 

\ e< l'- L 'q>/<p l 

(3) Determination of Bp: Equation (ll) gives 

£ _ _m_ _9_ 

5 m z ~ 

In order to determine from equation (ll ) , enter figure 4 

with W e(1 ^i = 7.0 and M^p = 3-0 to obtain ; then enter fig- 

* 1 / rn \ 


ure 


Then 


0 to obtain 


5 with (Ng^p)^ = 7.1^ and M^p = 3- 
is given by 

fi .) - ( -S-\ - f-2l\ = 0.6852 - 0.6202 = 0.0650 
Wp \ m lJ p \Pl /p 


( 


£P\ 

9 ih 


The boundary-layer thickness 6^ is then equal to 


Bp = = 1.272 X 10 meter 


Determination of Bp*: Equation (l4) gives 


81 = 1 jl 

S m 7 


Utilizing the information presented in step 3 yields 

{f\ ' 1 - ($) * 1 - °- S8 5 2 - °-5 148 
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and the value of 5 *^ is given by 


6*1 . Bj. 



= 4.002 X 10 “5 meter 


( 4 ) Change in thickness parameters 8 , 8 *, and 9 : 

(a) Change in 8 : Enter figures 19 (b) or 23 (b) (curves for 

^eq .,1 = 7 * 0 ^ with = 3*0 and a,j = 8 ° to obtain 


&d 


— = 0.716 or ^ = 0.968. 

81 Op 

(b) Change in 9 : Enter figure 20 (b) (curves for N e q^i = 7 *o) 

with = 3*0 and ctj = 8° to obtain ^ = I.O76. The 

momentum thickness change from station 1 to station E is 
assumed to be identical to the momentum thickness change from 
station 1 to station d. 

(c) Change in 8*; Enter figures 2 l(b) or 24 (b) (curves for 

N e q^i = 7*o) with Mj^i = 3-0 and aj. = 8° to obtain 

8* 8* 

= 0.866 or = 0.712. 

5 *1 5 eq 

(d) Change in N : Enter figure 17(b) (curves for N e g^i = 7 *o) 

Pj d / 

with M7 -> =3*0 and — - — = 3*008 (static-pressure rise 

v l,l V 

through incident -reflecting shock with an incident shock- 
turning angle of 8°) to obtain N d = 3. 34 . 

(5) The values of the thickness parameters 8 d , 9 d , and S* d at 

station d are: 

8 d = Si = 9.11 X 10 "^ meter 
a x §1 


9 - 


9 d = ®1 = 8.90 X 10 meter 


5 * d = S* n ^ = 3*468 X 10“5 meter 
a 1 8 *! 
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The values of the thickness parameters 
station E are: 


5 = S, = 1.233 X lO" 14 - meter 

cU J- 


eq 


5 


9 e q = 8.90 X 10“^ meter 


6* e g = 5*1 = 2.85 X 10“5 meter 

1 


* 


eq " - -l 5 * 
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(a) Boundary-layer shock patterns and pressure distributions. 
Figure 1.- Boundary- layer shock configuration. 
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(c) Analytical boundary- layer shock model. 
Figure L- Concluded. 
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(a) Mj between 1.25 and 5.75. 


ffl 


as a function of local free-stream Mach number M[ and profile index N. 
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(a) Shock waves. 


Figure 8.- Boundary-layer shock model for Yj = 61. 
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(a) Shockwaves. 

Figure 9.- Boundary-layer shock model for Yi > fij. 
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Figure 18.- Calculated boundary-layer thickness changes as compared with experimental data. 
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Figure 18.- Concluded. 





(a) Nj = 5 or R xe = 3.3 x l()5. 

Figure 19.- Curves for the boundary-layer thickness change between stations 1 and d. 






(c) Ni = 9 or R xe = 3.3 x 10?. 
Figure 19.- Concluded. 
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(c) Ni = 9 or 


R xe = 3.3 x 10?. 


Figure 21.- Concluded. 



















(a) Ni = 5 or R xe = 3.3 x io5. 

Figure 23.- Computed curves for the boundary- layer thickness change from upstream of the shock, station 1, 
to the equilibrium velocity profile station, station E, downstream of station d. T^t = 310° K. 




















(b) Nj = 7 or R xe = 3.3 x 10 6 . 
Figure 24.- Continued. 










(b) Variation of incompressible average friction coefficient with R xe - 
Figure 25.- I ncompressible average friction coefficient for turbulent flow. 
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